EIGENVALUES OF GUE MINORS 



KURT JOHANSSON AND ERIC NORDENSTAM 



Abstract. Consider an infinite random matrix H = {hij)o^ij picked from the Gaussian 
Unitary Ensemble (GUE). Denote its main minors by Hi = (/irs)i<r,s<i a^nd let the j:th 
largest eigenvalue of Hi be /i* . We show that the configuration of all these eigenvalues (i, fij ) 
form a determinantal point process on N X R. 

Furthermore we show that this process can be obtained as the scaling limit in random 
tilings of the Aztec diamond close to the boundary. We also discuss the corresponding limit 
for random lozenge tilings of a hexagon. 



This version of this article differs from the one pubUshed in Electronic Journal of Probability 
in that the errors listed in the separate erratum have been corrected. 

1. Introduction 

The distribution of eigenvalues induced by some measure on matrices has been the study 
of random matrix theory for decades. These distributions have been found to be universal in 
the sense that they turn up in various unrelated problems, some of which do not contain a 
matrix in any obvious way, or contain a matrix that does not look like a random matrix. In this 
article, we propose to study the eigenvalues of the minors of a random matrix, and argue that 
this distribution also is universal in some sense by showing that it is the scaling limit of three 
apparently unrelated discrete models. 

The largest eigenvalues of minors of GUE- matrices have been studied in |Bar01j , connecting 
these to a certain queueing model. It is a special case of the very general class of models analysed 
in |Joh03| . The large N limit of this model will yield the distribution of all the eigenvalues of a 
GUE-matrix and its minors. 

This process will turn out also to be the scaling limit of a point process related to random 
tilings of the Aztec diamond, studied in [JohOSaj and of a process related to random lozenge- 
tilings of a hexagon, studied in |Joh05b| . 

1.1. Eigenvalues of the GUE. Consider the following point process on A = N x M. There 
is a point at {n, fi) iff the n:th main minor of H, i.e. i?„, has an eigenvalue fj,. We will call 
this process the GUE minor process. A central result in this article is that this process is a 
determinantal point process with a certain kernel iiT*^^^. 

For details of what it means for a point process to be determinantal, see section [51 An explicit 
expression for this kernel is given in the next definition. 

Definition 1.2. The GUE minor kernel is 



-1 
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(r + j)! 
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where ^; s,ri) — when r < s and 



(r — s — Ij! 

-^^^ — y , ^ / it^^)-^-^-'e-' dt 

^ VF+Jy-i-s-j-iy-Jv 

for r > s. 

Here, hk{x) — 2~''^^{k\)~^^^T:~^^'^Hk{x) are the Hermite polynomials normalised so that 



hi{x)hj{x)e ^ dx = 5ij, 
hk = when k < and H is the Heaviside function defined by 

(l.f) Hit) = 



1 for i > 
for i < 0. 



Theorem 1.3. The GUE minor process is determinantal with kernel K'^'^^ . 
This will be proved in section |31 



1.4. Aztec Diamond. The Aztec diamond of size N is the largest region of the plane that is 
the union of squares with corners in lattice points and is contained in the region |a;| + |?/| < iV+1, 
see figure [TJ 

It can be covered with 2x1 dominoes in 2^(^+1)/^ ways, |EKLP92a[ lEKLP92b) . Probability 
distributions on the set of all these possible tilings have been studied in several references, for 
example |Joh05a[ IPro03| . Typical samples are characterized by having so called frozen regions 
in the north, south, east and west, regions where the tiles are layed out like brickwork. In the 
middle there is a disordered region, the so called temperate region. It is for example known that 
for large iV, the shape of the temperate region tends to a circle, see |JPS98) for precise statement. 

The key to analyzing this model is to colour all squares black or white in a checkerboard 
fashion. Let us chose colour white for the left square on the top row. A horizontal tile is of type 
N, or north, whenever its left square is white. All other horizontal tiles are of type S, or south. 
Likewise, a vertical domino is of type W, or west, precisely if its top square is white. Other 
vertical dominoes are of type E. 

In figure m tiles of type N and E have been shaded. Notice that along the line i = 1, there is 
precisely one white tile, and its position is a stochastic variable that we denote x\. Along the line 
1 = 2 there are precisely two white tiles, at positions x\ and a;| respectively, etc. In general, let x], 
denote the j-coordinate of the fc:th white tile along line i. These white points can be considered a 
particle configuration, and this particle configuration uniquely determines the tiling. It is shown 
in |Joh05 a that this process is a determinantal point process on = {1,2,..., iV}^, and the 
kernel is computed. 

We show that this particle process, properly rescaled, converges weakly to the distribution for 
eigenvalues of GUE described above. More precisely we have the following theorem that will be 
proved in sectional 

Theorem 1.5. Let /i* he the eigenvalues of a GUE matrix and its minors. For each N, let {x^j} 
be the position of the particles, as defined above, in a random tiling of the Aztec Diamond of size 




Figure 1. An Aztec Diamond of size 20 with N and E type dominoes shaded. 
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Figure 2. Lozenge tiling of a hexagon. 



N. Then for each continuous function of compact support : N x 



with < < 1, 



E 



= lim E 

N^oo 



N/2 



VW2 



)) 



1.6. Rhombus Tilings. Consider an (a, 6, c)-hexagon, i.e. a hexagon with side lengths a, b, c, 
a, b, c. It can be covered by rhombus-shaped tiles with angles tt/3 and 2tt/3 and side length 1, 
so called lozenges. The number of possible such tilings is 

a b 

nnn 

1=1 j= 



=1 fc=i 



i+j+k-l 
i + j + k-2' 



This formula was proved by Percy MacMahon (1854-1928), see |Sta99[ page 401] for historical 
remarks. 

Thus, we can chose a tiling randomly, each possible tiling assigned equal probability. A typical 
such tiling is shown in figure[2J Just like in the case of the Aztec diamond, there are frozen regions 
in the corners of the shape and a disordered region in the middle. It has been shown, that when 
a = b = c = N^ oo, this so called temperate region, tends to a circle, see jCLP98j for precise 
statement and other similar results. 

Equivalently, consider a simple, symmetric, random walks, started at positions (0,2j), 1 < 
j < a. At each step in discrete time, each walker moves up or down, with equal probability. 
They are conditioned never to intersect and to end at positions {c + b, c — b + 2j). Figure [31 the 
red lines illustrate such a family of walkers, and shows the correspondence between this process 
and tilings of the hexagon. These red dots in the figure define a point process. |Joh05b) shows 
that uniform measure on tilings of the hexagon (or equivalently, uniform measure on the possible 
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Figure 3. Tiled hexagon with sides a — 8, b — 5 and c — 10. The so called 
horizontal rhombuses are marked with a blue dot. 
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configurations of simple, symmetric, random walks) induces a measure on this point process that 
is determinantal, and computes the kernel. 

We will show, in theorem 15. 41 that the complement of this process, the blue dots in the figure, 
is also a determinantal process and compute its kernel. 

Let us introduce some notation. Observe that along the line t ~ 1, there is exactly one blue 
dot. Let its position be a;}. Along line t = 2 there are two blue dots, at positions and 
X2 respectively, and so on. All these are stochastic variables, and they are of course not 
independent of each other. 

We expect that the scaling limit of the process {x'j}ij, as the size of the hexagon tends to 
infinity, is the GUE minor process with kernel K'^^^. More precisely, let /i* be the eigenvalues of a 
GUE matrix and its minors. Then for each continuous function of compact support (j) : NxR — > M, 
< < 1, 



(1.2) E 



lim E 



nx', - N/2 
(l-0(z, ' , 1 )) 



We will outline a proof of this result by going to the limit in the formula for the correlation 
kernel, which involves the Hahn polynomials. A complete proof requires some further estimates 
of these polynomials. 

The GUE minor process has also been obtained as a limit at "turning points" in a 3D partition 
model by Okounkov and Reshetikhin jOR06j . We expect that the GUE minor process should be 
the universal limit in random tilings where the disordered region touches the boundary. 

Acknowledgement: We thank A. Okounkov for helpful comments and for sending the preprint [OR06) . 

2. Determinantal point processes 

Let A be a complete separable metric space with some reference measure A. For example R 
with Lebesgue measure or N with counting measure. Let Af(A) be the space of integer-valued 
and locally finite measures on A. A point process a; is a probability measure on M(A). For 
example, let x he a. point process. A realisation x{lo) is an element of A/(A). It will assign 
positive measure to certain points, {a;i(aj)}i<j<^(jj), sometimes called particles, or just points in 
the process. In the processes that we will study the number of particles in a compact set will 
have a uniform upper bound. 

Many point processes can be specified by giving their correlation functions, pn : A" — > E, 
n = 1, . . . , 00. We will not go into the precise definition of these or when a process is uniquely 
determined by its correlation functions. For that we refer to any or all of the following references: 
IDVJ881 Ch. 9.1, A2.1], |Joh05d[S5i00] . 

Suffice it to say that correlation functions have the following useful property. For any bounded 
measurable function (p with bounded support B, satisfying 

(2.1) y^Mk f p„(xi,...,a;„)d"x<oo 
the following holds: 

oo „ 

(2.2) E[[](l + 0(X,(W)))] = I] - / H^l) ■ ■ ■ (p{Xn)Pn{xi, Xn) d" X. 

Correlation functions are thus useful in computing various expectations. For example, if A is 
some set and xa is the characteristic function of that set, then 1 — E[]^(l — XA{xi))] is the 
probability of at least one particle in the set A. If the correlation functions of a process exist 
and are known, this probability can then readily be computed with the above formula. 
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We will study point processes of a certain type, namely those whose correlation functions exist 
and are of the form 

p„(a;i, . . . ,x„) ^det[K{xi; Xj)]i<ij<„, 

i.e. the n:th correlation function is given as a n x n determinant where A" : — )• M is some, not 
necessarily smooth, measurable function. Such a process is called a determinantal point process 
and the function K is called the kernel of the point process. 

Let , x^, . . . , ,. . . be a sequence of point processes on A. Say that assigns positive 
measure to the points {a;f^(w)}i<i<jv"(u;)- Then we say that this sequence of point processes 
converges weakly to a point process x, written x^ — >■ x, iV — ^ oo, if for any continuous function 
4) of compact support, < (/> < 1, 

'n{uj) 

l[il-<Pix,{L0))) . 

i=l 

The next proposition gives sufficient conditions for weak convergence of a sequence of determi- 
nantal processes in terms of the kernels. 

Proposition 2.1. Let ,. . . be a sequence of determinantal point processes, and 

let x^ have correlation kernel satisfying 

(1) — >■ K , N ^ oo pointwise, for some function K , 

(2) the are uniformly hounded on compact sets in A^ and 

(3) For C compact, there exists some number n = n{C) such that 

<lel[K^ {xi,Xj)]i<i^j<„i = 

if m > n. 

Then there exists some determinantal point process x with correlation functions K such that 
— >■ X weakly. 

Proof. We start by showing that there exists such a determinantal point process x. In fSosOOj, 
the following necessary and sufficient conditions for the existence of a random point process with 
given correlation functions is given. 

(1) Symmetry. 

Pkix^ii),- ■ • ,a;^(fc)) = Pk{xi,...,Xk) 

(2) Positivity. For any finite set of measurable bounded functions (j)k : A*"' — ^ R, fc = 0, . . . , M, 
with compact support, such that 

M 

(2.4) <t>0+J2 HX^^,■■■,X^,)>0 

for all (xi, . . . , xm) € /^^ it holds that 

N 

(2.5) ^o + Y (j)k{xi,...,Xk)pk{xi,...,Xk) dxi...dxn>0- 

k=i 

The first condition is satisfied for all correlation functions coming from determinantal kernels 
because permuting the rows and the columns of a matrix with the same permutation leaves the 



JV"(a;) 



(2.3) 



lim E 



n (l-^(a:f(a.))) 
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determinant unchanged. For the positivity condition consider the kernels . They are kernels 
of determinantal processes so 

M 

(2.6) ^^0+^/ 4>k{.xi,...,xu)det[K^{x^,Xj)]i<i,j<kdxi...dxn>0. 

As — > oo, this converges to the same expression with K instead of by Lebesgue's bounded 
convergence theorem with assumption (2). Positivity of this expression for all N then implies 
positivity of the limit. 

So now we know that x exists. We need to show that x^ — >■ x Take some test function 
: A — >■ K with bounded support B. For this function we check the condition in (j2.1l) . The 
assumption (3) in this theorem implies that the sum is a finite one. Also, ||0||oo < 1- Assumption 
(2) is that the functions are uniformly bounded, so in particular they are bounded on S^, 
so pk is bounded on B'^ . The integral of a bounded function over a bounded set is finite, so this 
is the finite sum of finite real numbers, which is finite. 

Therefore, for each N, by ([221), 



(2.7) lim E 



N- 



n(i-0(^f(^))) 



(2.8) ^^'^J2—r Y[Hx^)det[K^{x,,x,)]l<,^,<nd''X{x). 

Condition (3) guarantees that the sum is finite. Lebesgue's bounded convergence theorem applies 
because the support of 4> is compact and the correlation functions are bounded on compact sets. 
Thus the limit exists and is 

(2.9) =12^-r Y[(l^ix^)det[K{x,,x,)l<,,,<,-,d^X{x) 



(2.10) =E 



Ylil ~ cj,{x,{u;))) 



This implies that indeed x^ — >■ x, weakly, as A^ — cxo. □ 

3. The GUE Minor Kernel 

3.1. Performance Table. Consider the following model. Let {w{i, j)}(^ij^Qi2^, be independent 
geometric random variables with parameter q^. I.e. there is one i.i.d. variable sitting at each 
integer lattice point in the first quadrant of the plane. Let 

(3.1) G(M,A^)=max V w{i,j) 

{i,j)eiT 

where the maximum is over all up/right paths from (1, 1) to (M, A^). The array [G{M, N)]m, Net's 
is called the performance table. 

Each such up/right path must pass through precisely one of (M — 1, A^) and (Af, A^ — 1), so 
it is true that G{M, N) = max(G(Af - 1, N),G{M, A^ - 1)) + w{M, N). 

It is known from jBarOl) . that {G{N, 1), G{N, 2), . . . , G{N, M)) for fixed M, properly rescaled, 
jointly tends to the distribution of (/i}, /i^, . . . , fj^f^) as A^ — >■ oo in the sense of weak convergence 
of probability measures. We will show that it is possible to define stochastic variables in terms 
of the values w that jointly converge weakly to the distribution of all the eigenvalues /i* of 
GUE-matrices. 
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3.2. Notation. We will use the following notation from [BarOlj . 

(1) Wm,n is set of M X iV integer matrices. 

(2) WM,N,k is set of M x N integer matrices whose entries sum up to fc. 

(3) Vm ~ where the components of each element x are indexed in the following 
way. 

x\ 

„M-1 M 

(4) Cgc C Vm is the subset such that x^j_i > xy_\ > x'j. 

(5) Cgc,n are the integer points of Cgc- 

(6) Let p : Cgc be the projection that picks out the last row of the triangular array, 
i.e. p{x) = {x^ , . . . , x^j). Likewise, let q : Cgc.n — >■ N*^, the projection that picks out 
the last row of an integer triangular array. 



M-1 



Xi 



3.3. RSK. Recall that a partition A of fc is a vector of integers (Ai, A2, . . • ), where Ai > A2 > . . . 
such that Ai = fc. It follows that only finitely many of the A^is are non-zero. 

A partition can be represented by a Young diagram^ drawn as a configuration of boxes aligned 
in rows. The i:th row of boxes is Ai boxes long. A semi- standard Young tableau (SSYT) is a 
filling of the boxes of a Young diagram with natural numbers, increasing from left to right in rows 
and strictly increasing from the top down in columns. The Robinsson-Schensted-Knuth algorithm 
(RSK algorithm) is an algorithm that bijectively maps Wm,n to pairs of semi-standard Young 
tableau. For details of this algorithm, see for example |Sag01[ ISta99j . 

Fix a matrix w e Wm,n- This matrix is mapped by RSK to a pair of SSYT, (P{w), Q{w)). 
The P tableau will contain elements of M {1,2,..., M} only. Construct a triangular array 



M-1 M-1 

-I ■ ■ ■ a-M-i 

„M-1 M 
Xl ■ ■ ■ Xj^ 



where is the coordinate of the rightmost box filled with a number at most i in the j-th row 
of the P(?i')-tableau. This is a map from Wm,n to Cgc,n- 



3.4. A Measure on Semistandard Young Tableau. Consider the following probability mea- 
sure on Wm,n- The elements in the matrix are i.i.d. geometric random variables with parameter 

. Recall that a variable X is geometrically distributed with parameter g^, written X e Ge{q^) 
if P[X = fc] = (1 — q^){q'^Y , fc > 0. The square here will save a lot of root signs later. Such a 
stochastic variable has expectation a = q^ / (1 — q^) and variance h = q^/ (1 — q^)^- 

Applying the RSK algorithm to this array induces a measure on SSYT:s, and by the corre- 
spondence above, a measure on Cgc.n- Call this measure Tr^^^jy. The following is shown in 
IJohOOj . 

Proposition 3.5. Let Wu,n contain i.i.d. Ge{q^) random variables in each position. The 
probability that the RSK correspondence, when applied to this matrix, will yield Young diagrams 
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of shape A = (Ai , . . . , Am) is 
(3-2) Pf^,N:=^^^jJ—U 



Ml fJ-^jliN-M + jy. 

l<'i<j<A/ i=l ^ * ^ 

where fc = |A| A^. 

In other words, the measure t^^^^ integrating out aU variables not on the last row, is 
Pq^^ N- This, together with the following result characterizes the measure tt^^^j ^ completely. 

3.6. Uniform lift. Proposition 3.2 in jBarOll states that the probability measure tt^,^^^ ^, con- 
ditioned on the last row of the triangular array being A, is uniform on the cone q^^{X) := {x G 
Cgc,n : q{x) = A}. 

In formulas, this can be formulated as follows. 

Proposition 3.7. For any bounded continuous function : M x Z — >■ M of compact support, 

„[n(i+'^(*'^}))]=Efz7Ay E n(i+'^(*'^}))|'°s^w- 

where L(X) is the number of integer points in q^^{X). 
The number of such integer points is given by 



L(A)=n 



■ . J 

%<3 



3.8. GUE Eigenvalue measure. It is well know, see for example |Meh91| . that the probability 
measure on the eigenvalues induced by GUE measure on M x M hermitian matrices is 

p^y^{x,,...,XM)=^ n i^^-^^f n exp(-A?) 

l<i<j<Af l<i<M 

for some constant Zm that we need not be concerned with here. 

3.9. Uniform lift of GUE measure. |Bar01] shows a result for eigenvalues of minors of GUE 
matrices that is similar to the above result for partitions. He shows that given the eigenvalues 
of the whole matrix A = (Ai > • • • > Am), the triangular array of eigenvalues of all the minors 
are uniformly distributed in p^^(A) :~ {x G Cqc '■ p{x) — A}. Again we can write this more 
formally. 

Proposition 3.10. For any bounded continuous function (/) : M x R — > K. o/ compact support, 
the measure ttm^^ satisfies 

where Vol(A) is the volume of the cone p^^(A). 
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This volume is given by 

voi(A)=.n^^- 

1 — I 

This situation is then very similar to the measure t^^^^j above, in the sense that, conditioned 
the last row, the rest of the variables is uniformly distributed in a certain cone. 

3.11. Scaling limit. We are now in a position to see the connection between the measures 

„RSK J „GUE 

%2 M.Af ^^"^ ■ 

Proposition 3.12. Lef a E[w(l, 1)] = - 9^) and b := Yar[w (1,1)] ^ q'^ / {1 - q'^)'^ . Then 

for any bounded continuous function of compact support 4>, 

— aN 

IE.ouB[n(l + A^}))] = lim E^RSK [[](1 + <j,{i, = ))]. 

Proof. Plug in the expression for the right hand side in proposition 13.71 and for the left hand 
side in proposition 13 . lOl Stirling's formula and the convergence of a Riemann sum to an integral 
proves the theorem. □ 

3.13. Polynuclear growth. The measure tt^.^^ is a version of the Schur process and is a 
determinantal process on M X N, by |OR03| . We will use the following resuh from |JohQ3| . 

Proposition 3.14. The process {x^j} with the measure described in \3.4\ is determinantal with 
kernel 

(33) (rxsv)- ^ fdz fdw z wV^^ (1 - qwY (z ~ g)^-^^ 

For r < s, the paths of integration for z and w are anticlockwise along circles centred at zero 
with radii such that q < \w\ < \z\ < 1/q. For the case r > s, integrate instead along circles such 
that q <\z\ < \w\ < 1/q. 

This follows immediately from proposition 3.12 and theorem 3.14 in | Joh03j . 

Having now introduced the PNG- kernel, we can state the following scaling limit result. 

Lemma 3.15. Let a = q^ /(I — q^) and b = q^ /(I — q^)^ as above. The following claims are true 
for M fixed. 

(1) For r, s < M, 

V2WVif ™'/(r, [aN + ^V2bN\ ; s, [aN + r,V2bN\ ) — > K^'^^ir, C; s, r,) 

uniformly on compact sets as N oo for a certain function g 0. 

(2) The expression 

ipil^V2bNK^''^{r, [aN + ^V2bN\ ; s, [aN + 7]V2bN\ ) 
g{s,T],N) 

is bounded uniformly for 1 < r, s < M and ^, rj in a compact set. 

The proof, given in sectionlHl is an asymptotic analysis of the integral in p.3p . Now everything 
is set up so we can prove the main result of this section. 

Proof of theorem \1.3[ According to proposition 13.121 

(3.4) Kfr^Yli^ + = J™ ^^^i^ [11(1 + ^7=^))]. 
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The point processes on the right hand side of this last expression are determinantal. Their kernels 
can be written 

K^(r, s, r,) ■■= ^P^^V2bNK^''^{r, aN + CV2bN; s, aN + r,V2bN). 

for some function g that cancels out in all determinants, and therfore does not affect the corre- 
lation functions. By lemma [5. 151 these satisfy all the assumptions of proposition [5111 Thus, 
the point processes that these define converge weakly to a point process with kernel K'^^^. This 
implies that the measure on the left hand side of equation (13. 4p . i.e. irf}^^ is determinantal with 
kernel K'^^^. The observation that M was arbitrary completes the proof. □ 

4. Aztec Diamond 

The point-process connected to the tilings of this shape, described in the introduction was 
thoroughly analyzed in Joh05a) . The following result is shown. 



Proposition 4.1. The process {x'j} described in section is determinantal on A = N x N, 
with kernel given by 

(4.1) KU2r,.,2s,y)- ' ^^^--^--^(1 + 1/-)-- 



(27ri)2 7 z J w z^(l - z)''(l + 1/z)^-'^ z-w 

and reference measure fj, which is counting measure on N. The paths of integration are as follows: 
For r < s, integrate w along a contour enclosing its pole at —1 anticlockwise, and z along a 
contour enclosing w and the pole at but not the pole at 1 anticlockwise. For r > s, switch the 
contours of z and w. 

Based on this integral formula we can prove the following scaling limit analogous to that in 
lemma 13.151 

Lemma 4.2. The following claims hold. 
(1) 

ip^^/N/2KU2r, [N/2 + ^^/N/2\ ; 2s, LA/2 + i^^/N/2\ ) ^ if ^"^(r, C; s, r;) 

uniformly on compact sets as A —> oo for an appropriate function g ^ 0. 
(2) The expression 

^P^, VNj2K^{2r, [N/2 + ^^nJ2\ ; 2s, LA/2 + ^^nJ2\ ) 

is uniformly bounded with respect to A for (r, ^, s, rj) contained in any compact set in 
N X R X N X E. 

The proof is based on a saddle point analysis that is presented it section [6] We can now set 
about proving the main result of this section. 

Proof of theorem ] 1.5i By proposition 14. 11 the form a determinantal process with kernel K^. 
The rescaled process (a:* — A/2)/a/A/2 has kernel 

(4.2) A^(r,es,ry) := ipil^yW/2KU2r, N/2 + ^^N/2; 2s, N/2 + fj^N/2). 

By lemma the kernels satisfy all the assumptions of proposition 12.11 So they converge 
to the process with kernel K'^^^. □ 
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5. The Hexagon 

Consider an (a,b,c)-hexagon, such as the one in figure[31 First we need some coordinate system 
to describe the position of the dots. Say that the a simple, symmetric, random walks start at 
t — and y = 0, 2, . . . , 2a — 2. In each unit of time, they move one unit up or down, and are 
conditioned to end up at ?/ = c — &, c — & + 2, . . . , c — 6 + 2a — 2 at time t — b + c and never to 
intersect. One realisation of this process is the red dots in figure |31 At time i, the only possible 
y-coordinates for the red dots are {at + 2fc}o<fc<'yt, where 

[-t 0<t<b 
t-2b b<t <b + c. 




jt = \o + a~L b <t <- c at 

c, 



Let Aa.b,c = {{t, at + 2k) : 0<t<b + c, 0<k< -ft} be the set of all the dots, red and blue. 

5.1. A determinantal kernel for the hexagon tiling process. We now need to define the 
normalised associated Hahn polynomials, q^j^\x). These orthogonal polynomials satisfy 

N 

(5.1) E - Km, 



x=0 

where the weight function is 

^ ^ ' xlix + ay.{N + p - xy.{N - xy: 

They can be computed as 

~{a,P)f x_ { — N — l3)n{~N)n j-, f-n,n-2N-a-P-l.-x t\ 
In.N - ~(aR) , 3-^^21 -N-0.-N ' '-^Ji 

where 

T(a,/3)\^ {a + (3 + N-l- n)N+i 



(a + /? + 2iV + 1 - 2?i)n!(/3 + - n)!(a + A^ - n)!(iV - n)! ' 
For convenience, let a^ = |c — r| and b^ = \b — r\. jJohOSb] shows the following. 

Proposition 5.2. The red dots form a determinantal point process on the space h.a,b,c with 
kernel 

Ka,b,c{^^ + 2a;; s, a^ + 2y) = -</>r,s(ar + 2a;, as + 2y) 




(a + s - 1 - n)!(a + 5 + c - r - 1 - n)! a , -.-b a , \ , / \ 

<ln,t{^)<ln":{y)^r{x)^s{y), 



(a + s — 1 — rt)!(a + 6 + c — s — 1 — n)! 
where (j)r,s{x,y) =0 if r > s and 

(l>r,s{x,y) 



y—x-\-s—r 
2 



otherwise. Furthermore, 



' {{br + xy{-/r + ar - xy.)-'^ 0<r<b 
UJr{x) — ^ {xl{'-fr + Or — 2;)!)^"'- b < r < c 

(a;!(7r - x)!)-! c<r<b- 
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and 

Uy\{ls-y)\)-^ 0<r<b 

^s{y)^ liibs + yy.hs-yy.)-^ h<r<c 

{ {{br + y)!(7r + ar - yy.)-^ c<r<b + c. 

It follows that the blue dots also form a determinantal point process. To compute its kernel 
we need the following lemma. 

Lemma 5.3. 

s — r 

s—r-\-2y-\-as—2x~ar 
2 



{a + r-l-ny.ia + b + c-s-l-ny.^'''^'- [y)^r[x)u;,[y) 

when s > r. 

Proof. This proof uses the results obtained in the proof of 15.21 in [JohOSbi equation 3.25]. Define 
convolution product as follows. For f ,g : 1? ^ "L^ define (/ * : — )■ Z by 

(/ * 9)i.x, y) := ^ f{x, z)g{z, y). 

Let (t){x, y) 6x,y+i + ^x^y^i- Also let 

<t'*^{x,y) := 5^^y 
(t)*^{x,y) := (l>{x,y) 



Set 



1 



/n 



{a-k){j~ky.{a-l-jy 
n\ (n — 2a — 6 — c + 1)^ 



kj (-a - 6 + l)fe(-a)fe 



and finally let 



ll){n, Z) := ^ /„,m ^ Cj^ra<t){'2j, z) 
m— j— 

0o,i(", y) := V'l", y) 

The dual orthogonality relation to (15. ip is precisely 

(5.2) 5^ gtv;'-)(x)gt:,;"-)(y)c..(x)^.(2/) = 4.,. 

n=0 

By equation (3.25), (3.30) and (3.32) of the above mentioned paper, 

0o,r(f^, Or + 22:) A{a,b,c,r,n)q^^i;^''\z)ibr{z), 

where 
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Inserting this into the orthogonahty relation in (|5.2[) gives 

Ulr(x) 

Convolving both sides of the above relation with 0* gives 



CI/ • 0- C- 7^- 7^ 1 

n=0 V : 7 I 1 ; 

which, when the left hand side is simplified, gives 

„_Q ^l^i c, r, nj 

Invoking equation (j5.ll) again to simplify the left hand side and explicitly calculating the right 
hand side gives 

It is easy to check that 



A{a, b, c, s,n) / (a + s — 1 — n)\{a + 6 + c — r— 1 — n)! 
A(a, b,c,r,n) V {a + r — 1 — n)l{a + b + c — s — 1 — n)\ ' 

which proves the lemma. □ 
We now need to introduce the normalized Hahn polynomials <ln N\^)- These satisfy 

N 

(5.3) Y.&^^)&^^)^N-'\^) - 5rr.,n, 

where 



(5.4) w^^-'^\t) = 



t\{N~ty. 

Theorem 5.4. The blue dots form a determinantal point process on the space Aa^b,c with kernel 



Ka,b,c{r-,x] s,y) = 



n— — C30 

w/ien s > r, and 
Ka.b,ci^,x; s,y) = 



n=0 

otherwise. 
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Proof. It is well known that the complement of a determinantal point processes on a finite set 
with kernel K is also determinantal with kernel K = I — K , i.e. K(x, y) = 5x,y — K(x, y). 

Applying this result to our problem, we consider 6x.ySr,s — b c('''^' s,y). We now separate 
two cases. When s > r see that 



K{r,x; s,y) 



^ V (a + ^-l-^)!(a + & + c-s-l-n)!^"'^'- ^'-^^ [y)ujr(x)uj,(y) 

n— V ^ ^ ^ ^ 

is a candidate for the kernel for the blue particles. By lemma 1^751 this simplifies to 
K{r,x; s,y) = 



y (a + r — 1 — n)\[a + b + c — s — 1 — ny. ' ' '° 



For s < r we just get 
K{r,x; s,y) = 



We now exploit a useful duality result from |Bor02| . It states that 



&i^)\/y^^N'\-) = {-ir<j^j^L%ix).Jw^^''\x). 

Insert this into the formulas above and define the new kernel 

Ka^b.c{r,x; s,y) := (-1)^^'' ^/oJs{y)uJr{x)-'^uJr{x)Cjs{y)~'^K{r, x; s,y). 

This kernel gives the same correlation functions as K, since the extra factors cancel out in the 
determinants. The new kernel can be written as 



Ka,bAr, x; s, y) = y^u;^^'"'^\x)u;^;""''^(y)x 



{a + r-1- n)\{a + b + c-s-l- n)!'^'^'-"'^'-^^^'^^^-"'^^^^^' 



when s > r, and 



Ka,b,c{^, x; s, y) = - \/w'^;''''-\x)w^'';'''''\y) x 



a-l 



V- ja + s-l- n)\{a + b + c - r - 1 - n)! jt^^a^) („\ 
2^\l {a + r-l-ny.ia + b + c-s-l-n)\ '^T'-".^'- ^^^97.-n,7. ^V' 



n=0 



otherwise. 

The change of variables j := a — 1 — n puts these expressions on a simpler form, thereby 
proving the theorem. □ 
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5.5. Asymptotics. Let < p < 1 be some real number. Let a = IpN , (5 — 7(1 — p)N ^ 

i^[pN + ^2p{l-p)N{\ + ^-'^)x\ . Then 

(5.5) ^2p(l - p)N{l + l-^)^w^:f{i)q'-:f{i) ^ i-irV^hrXx) 

uniformly on compact sets in a; as — >■ 00. For completeness we give the proof of this result in 
the appendix. 

We would like to apply this with p = \ and 7 = 2 to our kernel and letting a = 6 = c — > cxd, 
i.e. we would like to take the limit 

K{r,^;s,f]) = 

lim {-3Nr-'^3N/lK^^Ur, 1^/2 + ^V^N/ll ; 1^/2 + vV^\)- 

j\ —a— 00 

The factor (— 3A^)^~^ cancels out in all determinants and is thus of no import. For s > r we get 

(5.6) K{r,^;s,v)^ E \/ W^"^^'^"'^^' 

j = -oo V ^ 

and formally, if we ignore the fact that this turns into an infinite sum, for s < r we get 

(5.7) K{r, s,v)^- £ ^|^±||/.,+,(0/i.+,(^)e-(«^+''^)/^ 
This expression can be simplified with the following lemma 

Lemma 5.6. Let H be the Heaviside function defined by equation above. Then, 

(5.8) 7^^(^ - y)''-'H{x -y)^J2 \j^-^K-k{y)K{x)e-y 

n—k 

rit^yr^-^er^'dt 

pointwise for x ^ y. 

The proof is given in section [B) 

In view of this result, the infinite series in 15. 71 converges and the kernel K is exactly the GUE 
minor kernel K'^^^. The interpretation of this is the following. The distribution of the blue 
particles, properly rescaled, tends weakly to the distribution of the eigenvalues of GUE minors 
as the size of the diamond tends to infinity, equation (jl.2p . The only thing needed to make this 
a theorem is appropriate estimates of the Hahn polynomials to control the convergence to the 
infinite sum. 

6. Proof of lemmas 

Proof of lemma \5.6l As the Hermite polynomials are orthogonal, there is an expansion of the 
function in the left hand side of (|5.8p of the form 



(6.1) 



00 

{x - yf-^H{x -y)^Y. c«(2/)^fn(x), 

n=0 
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where Hn is the n:th Hermite polynomial, as defined in for example |KS98) . and where the 
coefficients are given by 

1 f°° 

(6.2) c„(y) = -— ^ / (x ~ yf-^H{x - y)H,,{x)e--" dx. 

It is known that Hn{x) — i/„_i(a;)) for n > 1. Integration by parts and limiting 

the integration interval according to the Heaviside function gives 

{x-yf-^H.n{x)dx^ / {k~l){x-yf-^Hn-i{x)dx 
Jy 

For n > k, repeat this process k — 1 times to get 
For < n < fc, stop doing partial integrations when Hq is reached, giving 

Inserting (I6.3p and (I6.4p in (|6.ip and changing to normalized Hermite polynomials proves the 
lemma. □ 

Proof of lemma V3.15\ Assume first that r < s. By proposition 13.141 we have to consider the 
integral 

(27ri)2 J^^^ J^^^ w z-w ^,,V2WV (1 _ qzY {z - g)^^ ' 

where jr is a circle around the origin with radius r oriented anticlockwise, q < ri < r2 < l/q, 
and 

(6.5) fiz) - log(z - q) - (1 - g^)-! log z. 



(Here we have ignored the difference between aN + ^\/2bN and its integer part.) Note that 
f'{z) = gives z = l/q. This leads us to choose 

1 2 



ri 



and to deform 7^.2 to a circle F oriented clockwise around l/q with radius l/a^jN/2. The specific 
choice of radii are convenient for the computations below. Choose 



9 

Then, 



(r,e,A^) = 2-'-/2e-«'/\-«^ 



a^/N/2 



(6.6) 4^^^^^^2WVif _^ ^y2WVj ; s, [aA^ + Ti^/2bN\ ] 



2bN 



q \ V2bN 
ay^Ij/2 J (27ri) 



n2 



w z-w ^TiV2bN {l-qzY {z-q)^ 
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Parameterize 7^ by w{t) = rie**-^", —tt/En <t < tt/Em, En — q/ay^N/2. We have 

w{0) -q 



Re{f{w{0)) - f{w{t))) - In 



w{t) - q 



^1 f 2nq{l- cos ENt) \ 
2 ^'Iv (^1-9)^ J 

< -^In {1 + q^ {1- cos ENt)) , 

for large enough. Since cosx < 1 + when \x\ < tt, the last expression is 

< In (1 + q^EltyS)) < -Ct^/N 
for \t\ < IT /En, where C > is a constant depending only on q. Hence, 

(6.7) ReiV(/(u;(0)) - f{w{t))) < ~Ct^ 

for \t\ < Tr/En with C > 0. 

In the right hand side of (|6.6p we make the change of variables 



(6.8) z = z{u) = l/q-u/ay/N/2 

with u on the unit circle oriented anticlockwise. We obtain the integral 



(6.9) \/2'^-'-e'''-«' 



du 



dt- 



1 



2iqVb 

{2mYa J^^ ""J-^/E^ a^/Nj2{z{u) - w{t)) 



^N{f{z{u))-f{w{t))) 



{qw{t)) 



r]V2bN 



M 



g \ (1 - qw{t)y {w{t) - q) 

" iqz{u))i^ \ a^/Nj2 J (1 - qziu)y iz{u) - q)^' 

Note that qVb/a = 1. Also, 

(6.10) f{l/q + h) = f{l/q) - + 0{h^) 
for \h\ small. Hence, for N sufficiently large, 

(6.11) Nifiz{u)) ~ f{w{0))) = -u^ + A + hNiu)/VN, 
where hNiu) is bounded for \u\ = 1. We have 



(6.12) 



{qw{t)) 



riV2bN 



riV2bN 



-fV2bN 



qu 



a^a^/jr/2^ 



By the inequality (1 + a:/n)" < for x > —n, n> 1, the right hand side in (|6.1ip has a bound 
independent of TV. We also have 



(6.13) 

for M € 71, |i| < tt/En, and 
(6.14) 



ay/N/2\z{u) ~ w{t)\ > 1 
\ (1 - qw{t)r {w{t) - q)^ 



a,/Nj2 {l-qz{t)Y {z{u)-q) 



M 



for u G 71, \t\ < tt/En, by (|6.8p and the definition of w{t). 
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It follows from (g^l), (ICTTI) . (|ST^ . (IBT^ and ((Ol) that the part of the integral in 
where the t-integration is restricted to N-'-/^ < \t\ < tt/En can be bounded by 



/|t|>JVl/3 

which goes to as iV — oo. When \t\ < N^^"^ we have 



(6.15) 



(1 - qw{t)) - (2 - it) 



< 



C 



iVl/6- 



Hence, for \t\ < N^^^ we have the bound 



(6.16) 



(1 - qw{t)r (wit) - q) 



I (1 - qz{u)y iz{u) - qY^ 



and we see that the part of the integral in (|6.9p where \t\ < N^^^ has a uniform bound for ^, rj 
in a compact set. This proves claim (2) in lemma [5. 151 for r < s. 

It also follows from (HJI), dUl), (I5T0)) . ((6ll|) . (EH]), ((05t . ((6l6)) and the dominated con- 
vergence theorem that the integral in (j6.9l) converges to 



(6.17) 



du / (i< 



1 



it — u 



u'^ (2-it f-2v(u-it) ~ ^^)^ 



Now let w = 2 — ii. Then we should integrate v along the line Re v = 2 from minus to plus 
infinity, call this contour F'. We obtain the integral 



du 



dv ■ 



2(7ri)2 

Expand {v — u)~^ as a geometric series. This turns the expression into 

2^u-u^ 



(6.18) 



2(7ri) 



n2 



and we recognize the classical integral representations of the Hermite polynomials. The expres- 
sion now becomes 



(6.19) 



fe=0 



'(s-fc- 1) 



(r — fc - 



which proves claim (1) in the lemma in the case r < s. 

We now turn our attention to the case r > s. Deforming the ly-contour through the z-contour 
in (|3.3p . we get the same integral as above save for a residue that we pick up at 2 = w. This is 



1 



We see that the argument above goes through for the remaining integral also when r > s 
until (|6.18p . From there the terms fc = 0, . . . , s -f 1 give (j6.19l) as before. In terms k ~ s + 2, 
. . . , r -|- 1 we instead evaluate the ii-integral using the formula 



(6.21) 



1 



pV ~2rjv 



■ dv 



2" 



TTi Jy f" -\/7r(n — 1)! 7^ 

which is valid for n > 1 and will be proved below. That accounts for the terms j ~ — r, 
in definition 11.21 



s + 1 
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Using the well known formula 

1 fn + k — 1 

(1 - xY ~ ^ 

^ ' fc=0 



the integral in (|6.20l) becomes 
It is readily solved as 



if y > a;. 



^\q^-y iiy<x 

x~y ^ — 



With our rescaling, x = aN + £^\/2bN and y — aN + rj\/2bN and the factors g{r, N) / g{s,r], N) 
we see that the integral in (I6.20p is 



(6.22) V2p!^^n^q(^~i)^y^{ — ^=)'-''V2WVx 

a-\/ N/2 



r{{^ - f])V2bN){r -s - 1)! 
where H is the Heaviside function. As A'^ — ;> oo we get the limit 



(6.23) ^g^2_j2 2r-s(|_^ H{^-r]), 

(r — s — 1)! 

at least for £, i]. The case ^ = ij is a set of measure zero and is not important. Together with 
the result for the double integral this completes the proof of claim (1). It remains to show the 
estimate in claim (2) in this case. But this is easy. The expression in (|6.22p is the exact solution 
of integral (j6.20p . and since this is bounded in N for ^, 77 in a compact set, claim (2) follows. 
It remains to show the formula (j6.21l) . Observe that, by repeated partial integration, 

,,_ (n-l)!e-^''- 



r—s— 1 



(6.24) (^_^)"-ie-^^-de 



(2t;)" 

if 77 > and w G F'. So in this case the left hand side of our formula can be written 
(6.25) 



ni{n 



Here we can change the order of integration and evaluating the Gaussian integral gives the right 
hand side of (l63T|) . 

When ?7 < 0, make a change of variables v ^ —v. The left hand side of (|6.2ip becomes 
(6-26) ^ / ——dv, 



where F" is parameterised v — —2 + it, t — —00 t-^ 00. By deforming the contour F" into F' we 
get 

(6.27) / '---dv+^-S- / '-——dv, 
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where 7 is a circle around the origin. The right term can be evaluated using the the result for 
rj > above. The left term can be rewritten using the equality between the two integral formulas 
for the Hermite polynomials mentioned above and we obtain 



(6.28) 



2" 



2" 



V^(n- 1)! J_ 



V^(^-i)!7-c 

which proves formula (I6.2ip for 77 < 0. □ 
Proof of lemma \J7^ Assume first that r < s. By proposition UTT] we have to consider the integral 

dw 1 



dz 



^NU(z)-f{w))'_ 



w z ~ w 



^r+y/N/2i (1 - zY (f + wY ' 



where 7^ is a circle around —1 with radius r oriented anticlockwise, 1 < ri < r2 < 2 and 

/(z) = ilnz-ln(f + z). 



(Here we have ignored the difference between N/2 + and its integer part.) In the proof 

of lemma [5.151 we could chose the contours of integration as circles centred at the origin. This 
cannot be done here. 

Note that f'{z) — gives z = 1. This leads us to choose 

2 

n = 2 ^= 

and to deform 7^2 to a circle T oriented clockwise around 1 with radius 1/ yjN/S. The specific 
choice of radii are convenient for the computations below. Choose 



Then, 



(6.29) ^p^JWj2K\r,[N/2 + £,^/Nj2\- s,YN /2 + ^^/nJ2\) 



Parameterize 7^ by 
(6.30) 



(27ri)2 



dz 



dw 1 



w z — w 



w{t) = — 1 + rie 



itEr. 



for -tt/En < t < tt/En, En ^ \I^J~Nl2. We have 



Re(/(«;(0)) - /(«;(t))) 



■In 



u;(0) 



,(t) 



--lull 



1 



2ri(l — cosE^t) 



in - ly 



1 



< --In 1 + -(1 - cosENt) 



for large enough N. Again cosx < 1 — when |a;| < tt, the last expression is 



< -- In (1 + Ej^t^/m) < -Ct^/N 
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for \t\ < n/Epf, where C > is an absolute constant. Hence, 

(6.31) Re{f{w{0)) - f{wm < ^Cf/N 

for \t\ < tt/En, with C > 0. 

In the right hand side of (|6.29l) we make the change of variables 

(6.32) z = z(u) = l-u/VWs 

with u on the unit circle oriented anticlockwise, denoted 7. We obtain the integral 



(6.33) V7V«-'-e')'-« 



27r2 



du 



n/E„ 



dt- 



1 



,N{f(z{u))-f{w{t})) 



7 J-TT/Biv Z{u)-W{t) 

{w{t)Y-^+'^V^ (1 - wit)y (1 + z{u)y 



(z(m))'-+«v^ {1 - z{u)y {1 + w{t)y 

Note that 

(6.34) f{l + h) = f{l)-hy8 + 0{h^) 
for small \h\. Hence, for N sufhciently large 

(6.35) Nif{z{u)) - fiw{0))) ^~u^ + A + hNiu)/VN, 
where h^lu) is bounded for \u\ = 1. We have 

(6.36) 



for some constant C > depending on r. We also have 
(6.37) 

for M S 7 and \t\ < tt/En, and 



yA78|z(7i) -«;(<)! > 1 



(6.38) 



^^^ {i-w{t)y {i + z{u)Y 
{1 - z{uy- {I + w{t)y-^ 



for M e 7, |i| < tt/Stv, by (jO^ and (lOO)) . 

It follows from (lOTl) . (jOSl) . (jOS]) . (|071) and (jOg)) that the part of the integral 
where the t-integration is restricted to N'^/'^ < \t\ < tt/En can be bounded by 



'\t\>N^/^ 

which tends to as TV —> 00. When \t\ < N^^"^ and u € 7, we have 



(6.39) 

and 
(6.40) 



< c 



(z(u))'^+«V^/2 

{i-w{ty'' {i + ziuy 



< c. 



(1 - z{uy^ {i + w{ty'-^ 

where C depends on s, r and 77 but is independent of N. 

Hence, we see that the part of the integral in (j6.33p where \t\ < N^^"^ has a uniform bound 
for ^ and 77 in a compact set. This proves claim (2) for r < s. 



24 



KURT JOHANSSON AND ERIC NORDENSTAM 



It also follows from (|OT|) . (lO^ . (jOS)) . (jOTl) . (jOg)) . dOg]) . (|On)) and the dominated 

convergence theorem that the integral in (|6.33l) converges to 

^2s-rei'-i-^ f du [ dt- \ ^(2-^tf-2{2-^t)r,^-u-+2u^ [2 - itf 

2n^J^ {2-it)-u u- 

which is exactly the integral in (|6.17p . This proves claim (1) in the lemma in the case r < s. 

For r > s we can deform the contours one through the other to get the same integral as we 
solved above. On the way we pick up the residue of a pole at z = w. It is 

(6.41) ^ / ( 

2ni L, w \1 — w 



where x = lN/2 + ^a/ N/2\ and y = lN/2 + ■qy'N/2\ . The argument above goes through for 
the remaining integral also when r > s. We see that if 77 > ^, then x — y —oo and this last 
integral is zero. For simplicity, let fc = r — s and (3 — {x — y)/y^N/2. The coefficient in front of 
in the expansion of [(1 + — w)]'^ is 

1 f dw _j f 1 + 'w\'' \^ f k\ fj + i — 1 



One then sees that the i = k term dominates when N is large. 

fc-i 



(6.42) 



i:(:)C"r')(-i) 



k—irxi 



< CN 



k-l 



i=0 

Keeping only the i = k term and plugging in our rescaling and the factors g(r,^, N)/g(s,ri, N), 
we see that the integral in (j6.41|) is 



(6.43) Vev'-i'N^-r^/N/22^-^(^'^^^^'^^'^'^) 

\ r — s — 1 / 



for ^ > 77. When ^ rj this tends to 



(r — s — 1)! 

as iV — > 00 which together with the corresponding result for the double integral settles claim (1) 
in the case r > s. 

Claim (2) in this case follows from the corresponding result for the double integral. (|6.42l) and 
the boundedness of the expression in (|6.43p . □ 

Appendix A. Asymptotics for Hahn polynomials 
The Hahn polynomials, as they are defined in |KS98| . satisfy 

E ^ " Qm{x- a, /3, A^)Q„(x; a, /3, N) = {d^:;S^)H,,m 

where 

(,(c.,P)^2 ^ (-!)"(» + g + /3 + l)Af+l(/3 + l)nn! 

^ "'^ ^ {2n + a + P + l)[a + l)„(-iV)„Af! ' 
The Hermite polynomials are defined as usual: 

(A.l) ^ / Hn{x)H„,{x) dx - T'n\5nm. 

With this notation, the following well known limit theorem holds. 



r—s— 1 
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Theorem A.l. Let < p < 1 and 7 > 0. Let x ^ [pN + a:^2p(l - p)N{l + 7-1)] and 



fn,N = i-iy\l[^)2-nl^ P ^ ( ^ 



n J \1 — p J \l + 7 

(A.2) Er,{x) = /„,jvQ„(£; WN. 7(1 - P)N, N). 

Then 

uniformly on compact sets. 

Proof. The idea is induction on n. To start with, Qq^tj, a, jS, N) ~ 1 and we actually have 
Eo{x) = Hoix). For n = 1, 

Qi{y,a,P,N) = 1 - 7-^7^77^2: 
[a + IjiV 

so 



= ---=Hi{x) + 0{^/n^). 

Now assume that the theorem is true for n and n — 1 . We wish to show that it is true for n + 1 . 
There are three term recursion formulas for both Hahn and Hermite polynomials. Let 

(n + a + + l)(n + a + 1){N - n) 



Ar 

Cn = 



{2n + a + 13 + l)(2n + a + /3 + 2) 
n(n + a + /? + iV + l)(n + /3) 



(271 + a + /^)(2n + a + ^ + 1) ' 
Then 

(A. 3) AnQn+l{x) = {An + C„ - x)Qn{x) ~ C„Q„_i(x) 

(A.4) Hn+i{x) = 2xiJ„(a;) - 2niJ„_i(a;). 

Solving (|A.2p for Qn and inserting into (jA.3p gives after some simplification the following: 

(A.5) En+l{x) = — 1 + 1- - 1- ^ri{x) - —En-.l{x). 

Jn,N \ An An J Jn-1,N An 

Observe that under our scaling, 

An^pN + OiN-^) 

7 

Inserting this into equation (jA.Sp and doing some manipulations gives 

En+i{x) = {2x + 0(A^"i/2)^ _^ (2n + 0{N-'^^)^ En^i{x), 

which with our induction assumption is 

= 2xHn{x) + 2nHn-i{x) + OiN-^'"^) 
^Hn+i{x)+0{N-^/^). 
This completes the proof. □ 
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Applying Stirling's approximation to d^'^, fn,N and the weight function w^^'^\x), it is easy 
to show that 

Corollary A.2. ^5 before, x = pN + x^2-p{\ - p)N{l + 7-1). 

(A.6) ^2p(l - p)Nil + l-')ql:fiS:)^Jw"/{x) (-l)"/.„(:c)e--V2 

as TV -> 00 ifa/N pj and f3/N -> (1 - p)^. 
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